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1.  Introduction 


The  problem  of  bonding  two  surfaces  is  important  in  many  aspects  of  Army  technology.  The 
need  for  bonding  may  arise  when  two  dissimilar  materials  require  mutual  attachment,  or  when 
geometric  complexity  and/or  cost  preclude  the  manufacture  of  an  item  as  a  single  piece.  Bonding 
requirements  may  call  for  a  permanent  adhesive  attachment  or,  alternately,  they  may  require  the 
ability  to  debond  the  items  on  demand,  as  in  the  case  of  tape  or  velcro.1  In  some  cases,  the 
substrate  materials  being  bonded  may  be  stiff,  relative  to  the  bonding  material,  while  in  other 
cases,  the  substrate  materials  may  be  wholly  flexible  ( e.g .,  tape,  velcro).  This  report  focuses  on 
bonds  involving  flexible  substrates. 

A  successful  program  in  adhesive  design  requires  a  knowledge  not  only  of  how  the  intact  bond 
successfully  performs,  but  also  of  how  and  when  the  bond  will  fail.  This  report  is  one  in  a  series 
of  reports  (/,  2)  in  which  we  examine  the  nature  of  bonding  between  two  surfaces,  and  the  energy 
interchange  involved  with  the  creation  of  new  free  surfaces  (i.e.,  fracture).  In  an  effort  to  reduce 
the  problem  to  its  essentials,  we  attempt  to  formulate  scenarios  that  are  maximally  simplified.  In 
this  way,  the  problem  can  be  more  fully  understood  without  being  lost  in  mathematical 
complexity.  In  this  report,  a  notional  experiment  is  described  and  analyzed,  involving  the 
progressive  fracture  of  a  one-dimensional  (1-D)  idealized  bond  under  controlled  loading 
conditions. 


2.  The  Idealized  Problem 


In  figure  1,  we  display  the  closed  system  that  we  intend  to  address  in  this  report.  It  consists  of 
two  cables  that  have  been  bonded  together,  which  are  being  separated  in  a  controlled  manner. 

The  word  “cable”  is  used  here  to  convey  one-dimensionality,  though  other  arrangements  that 
physically  approximate  one-dimensionality  would  apply,  such  as  velcro  straps  or  adhesive  tape. 
The  goal  of  such  a  notional  experiment  is  to  study  the  manner  and  degree  to  which  a  given  load 
mg  will  bring  about  the  fracture  of  the  bond,  altering  the  load  angle  90. 

1  Velcro  is  a  hook-and-loop  fastener  system,  which  is  a  registered  trademark  of  Velcro  Industries  B.  V.  In  this  report, 
we  use  the  lowercase  term  “velcro”  to  refer  to  generic  hook-and-loop  fasteners  in  the  style  of  Velcro. 
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Figure  1 .  The  system  under  consideration,  characterized  by  two 
bonded  cables  being  debonded  under  an  applied  load. 


In  figure  1,  point  P  is  a  rigid  attachment  point  for  the  bonded  cable  system.  Points  H  and  —H 
are  rigid,  frictionless  pins,  over  which  the  debonded  cables  respectively  pass.  The  load  applied  to 
the  cables,  brought  about  by  mass  m  under  the  influence  of  gravity  g,  can  and  will  be,  hereafter, 
characterized  by  a  magnitude  of  force  equal  to  T0. 

For  the  analysis,  we  focus  on  a  subset  of  figure  1,  as  shown  in  figure  2.  In  addition  to  the  zoom, 
figure  2  has  been  rotated,  for  convenience,  90°  clockwise  relative  to  figure  1 .  Here,  one  half  of  a 
symmetric  system  is  depicted  (/.<?.,  one  of  the  cables),  because  of  planar  symmetry  across  the 
x-axis.  The  idealized  cables  are  assumed  to  be  massless,  linearly  inextensible,  yet  fully  flexible 
( i.e .,  incapable  of  supporting  a  bending  moment).  It  is  important  to  remember  when  considering 
figure  2  that  the  cable  is  not  attached  to  point  II.  but  merely  directs  the  applied  load  T0  through 
point  H . 

The  artifice  of  planar  symmetry  was  introduced  to  this  problem  merely  to  ensure  that  we  were 
describing  a  physical  system  in  which  the  deforming  bonds  applied  force  purely  in  the 
^/-direction.  Mathematically,  and  with  no  loss  of  generality,  we  may  solve  the  half-problem 
actually  shown  in  figure  2  by  considering  the  cable  to  lie  upon  a  rigid  substrate,  under  the 
influence  of  a  perpendicular  force-field  of  specified  characteristics.  Note,  therefore,  the 
descriptive  similarity  to  the  thermomechanical  problem  described  in  reference  1. 
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Figure  2.  The  +y  configuration  of  an  idealized  experiment  [the  —y  configuration 

(not  shown)  is  symmetric  about  the  x-axis],  in  which  two  inextensible  but 
flexural  cables  that  have  been  bonded  and  attached  to  the  laboratory  at 
point  P  are  “unzipped”  by  applying  an  ever-increasing  tension  To  through 
points  H  and  —H  (not  shown). 


In  this  analysis,  we  assume  the  bond  provides  a  fixed  resistive  force  (/o)  per  unit  length  of  cable 
up  to  the  point  at  which  the  bond  breaks,  characterized  as  when  the  ^-displacement  exceeds  the 
critical  value  for  fracture,  u0.  Such  process-zone  mechanics  could  be  representative  of  several 
types  of  bonding  models,  for  example,  a  rigid-plastic  bond,  an  electric  field,  or  even  a 
gravitational  field.  Nonetheless,  for  a  bond  resulting  from  a  continuum  adhesive,  the  model  is 
highly  simplified,  in  that  transverse  (e.g.,  Poisson)  influences  on  the  stress  state  are  ignored. 

For  points  x  <  Sf,  the  cable  is  either  initially  unbonded  or  has  become  debonded  and  remains 
under  load  T0,  applied  at  an  angle  0o  with  respect  to  the  x-axis.  The  angle  60  will  change 
(i.e.,  decrease)  with  the  progression  of  the  bond  fracture,  as  the  load  level  T0  is  increased. 

For  cable  in  the  domain  Sf  <  x  <  Si,  the  deforming  bond  remains  intact,  but  under  resistive 
stress.  The  restoring  force  per  unit  cable-length,  given  by  /0,  is  independent  of  the 
^-displacement,  u.  The  tension  in  the  cable,  T(x),  will  not  remain  fixed  in  this  region,  since  the 
stressed  bond  is  applying  forces  to  the  cable  which  are  not  strictly  perpendicular  to  it.  This  part 
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of  the  problem  domain  is  of  £ -length  A,  though  the  pathlength  of  cable  in  this  domain  will  always 
exceed  A  (in  fact,  it  will  always  be  greater  than  or  equal  to  ^ A2  +  Uq  ). 

For  the  bonded  cable,  we  admit  the  possibility  that  beyond  a  certain  domain,  x  >  Si,  the  //-forces 
establish  a  precise  balance,  such  that  the  cable  displacement,  u,  remains  exactly  zero.  This 
domain  of  the  cable  is  still,  however,  subject  to  axial  tension  in  the  x-direction,  given  by  T~^. 

The  idealized  problem,  therefore,  consists  of  these  three  domains:  (1)  x  <  Sj-\  (2)  Sf  <  x  <  S, ; 
and  (3)  x  >  Si.  While  the  problem  envisions  that  the  stress  states  over  these  domains  and  indeed 
the  domain  limits  themselves  will  change  with  increasing  load,  T0,  the  problem  we  solve  here  is 
quasi-static,  in  which  the  load  T0  is  increased  slowly  and  incrementally,  allowing  progressive 
equilibrium  states  to  establish  themselves  throughout  the  system.  In  this  manner,  we  intend  to 
establish  a  relation  that  links  the  magnitude  of  the  applied  tension,  T0,  to  the  resultant  level  of 
debonding  as  characterized  by  domain  boundaries  S /  and  S',  (alternately  expressed  as  a  function 
of  60),  in  terms  of  the  relevant  material  parameters  (e.g.,  f0  and  u0). 


3.  Formulating  the  Boundary  Conditions  and  Governing  Equation 


3.1  The  Domain  x  <  S'/:  The  Detached  Zone 

The  domain  x  <  Sf  represents  that  portion  of  cable  that  was  not  bonded  to  begin  with,  as  well  as 
the  additional  portion  that  debonds  from  the  application  of  force  T0  upon  it.  The  only  force  acting 
upon  the  cable  in  this  domain  is  the  tension  applied  through  point  H.  Equilibrium  considerations 
dictate  that  the  cable  in  this  domain  of  the  problem  remain  straight  and  under  constant  load,  T0. 
There  is  no  requirement  that  the  cable  remain  forever  fixed  at  point  H.  Indeed,  quite  the  contrary 
is  required. . .  namely,  that  the  cable  be  pulled  through  point  H  as  the  load  T0  is  increased. 

3.2  The  Domain  x  >  S'*:  The  Bonded  Zone 

At  the  far  end  of  the  cable,  for  x  >  S',  out  to  point  P,  the  situation  is  likewise  simple  to  analyze. 
The  cable,  in  this  part  of  the  problem  domain,  lies  upon  the  x-axis.  The  force  applied  where  the 
cable  is  rigidly  attached  to  point  P  is  given  as  T^.  Force  equilibrium  dictates,  that  at  all  points  in 
this  part  of  the  problem  domain,  the  cable  tension  remain  fixed  at  T^. 

Furthermore,  by  considering  the  complete  length  of  cable  for  x  >  0  as  a  free  body,  and  knowing 
that  all  external  forces  upon  it  in  the  intermediate  region  (Sf  <  x  <  Si)  are  perpendicular  to  the  x 
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axis,  a  force  equilibrium  on  the  cable,  applied  in  the  ^-direction,  leads  one  to  conclude  that 


=  T0  cos  60 


(1) 


3.3  The  Domain  Sf  <  x  <  St :  The  Process  Zone 

We  now  proceed  to  analyze  the  most  important  part  of  the  problem — the  process  zone. 

3.3.1  The  Boundary  Conditions 

This  part  of  the  problem  domain  represents  the  key  part  of  the  problem:  the  deforming  bond  on 
its  way  toward  eventual  fracture.  We  may  conclude  several  salient  boundary  conditions  on  this 
domain.  Because  there  are  no  concentrated  forces  applied  at  St  and/or  Sf  (only  a  distributed 
force  along  the  entire  domain),  equilibrium  applied  across  the  domain  boundaries  will  require  a 
continuity  of  cable  tension,  displacement,  and  slope. 

The  cable  tension  (T)  changes  along  the  length  of  cable  in  this  part  of  the  problem  domain. 
However,  the  tension  must  be  T  =  T0  at  x  =  Sf  and  T  =  T^  at  x  =  St,  in  order  to  match  the 
adjacent  domains.  Since  the  coordinates  and  slopes  of  the  cable  must  also  match  the  adjacent 
domains,  this  requires  that,  at  x  =  Sf,  the  ^-displacement  u  =  u0  and  v!  —  —  tan  0O.  At  the  other 
end  of  the  domain,  x  =  St,  we  must  satisfy  the  requirement  that  u  —  0,  and  u'  =  0. 

3.3.2  The  Governing  Equation 

To  develop  the  governing  relations  for  the  process  zone,  we  must  satisfy  the  equilibrium 
requirements  for  an  arbitrary  section  of  cable  in  that  zone,  such  as  that  depicted  in  figure  3. 

In  the  ^-direction,  the  only  forces  arise  from  the  tension  in  the  cable,  and  therefore 

—  (Tcos0)  =  O  .  (2) 

ax 

From  this,  we  conclude  that  T  cos  6  must  remain  constant,  and  that  constant,  from  matching  the 
conditions  at  the  end  of  the  domain,  is  given  by 

T  cos  6  =  T^  =  T0  cos  0O  .  (3) 

In  the  y-direction,  the  situation  is  slightly  more  complex,  in  that  the  independent  variable  will  still 
be  taken  as  x,  even  as  the  forces  being  summed  are  y-forces.  If  F  represents  the  integral  along 
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Figure  3.  Force  balance  on  an  element 
of  bonded  cable  in  the 
idealized  process  zone. 


the  £ -direction  of  the  distributed  bond  force  /0,  then  the  condition  of  equilibrium  in  the 
t/-direction  becomes 


-^-(Tsin#  +  F)  —  0  .  (4) 

ax 

One  must  note,  however,  that  the  bonds  which  produce  the  distributed  force  /0  are  distributed 
along  the  cable  length  ds.  Therefore,  as  the  angle  0(x)  changes,  the  density  of  bonds  along  dx 
changes  as  well.  Thus,  the  quantity  dF/dx  must  be  considered  as  dF/ds  ■  ds/dx.  The  quantity 
dF/ds  is  merely  the  distributed  force  of  the  bond,  namely  /0,  while  ds/dx  is  obtained  from  the 
geometry  as  sec  8,  alternately  expressible  as  \/l  +  u'2.  Using  equation  3,  we  can  reexpress 
T  sin  8  as  tan  8,  alternately  expressible  as  —T^u'. 

Combining  all  these  substitutions  into  equation  4  yields  the  following  governing  relation: 

u  -  AVl  +  u'2  =  0  .  (5) 

^  OO 

Equation  5  should  be  immediately  recognized  as  being  associated  with  the  so-called  “catenary” 
solution,  which  is  the  curve  that  governs  the  geometry  of  a  free-hanging  chain  or  cable  under  the 
force  of  gravity  (indeed,  the  word  catenary  comes  from  the  word  catena ,  Latin  for  chain).  Such  a 
connection  should  not  be  surprising  when  one  considers  that  the  distributed  force  of  the 
deforming  bonds,  /0,  constant  per  unit  length  of  cable,  is  directly  analogous  to  the  distributed 
force  of  a  cable’s  weight  under  the  force  of  gravity. 
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It  may,  therefore,  be  instructive  to  visualize  the  solution  to  this  mechanics  problem  in  terms  of  the 
analogous  gravitational  problem,  in  which  a  long,  heavy  chain  is  laid  in  a  straight  line  upon  a 
table  (the  rigid  substrate).  As  a  free  end  of  the  chain  is  lifted  to  a  given  height  above  the  table 
(analogous  to  u0 ),  the  shape  in  which  the  hanging  chain  establishes  itself  will  vary  with  the  level 
of  tension  applied  (itself  being  a  function  of  the  vertical  (lifting)  and  lateral  (stretching) 
components  of  that  force).  The  derivation  of  the  catenary  solution  is  provided  in  the  appendix, 
for  the  reader’s  convenience. 


4.  The  Solution 


4.1  The  Catenary  Solution  Applied  to  the  Domain  Sf  <  x  <  Si 

To  solve  equation  5,  let  us  first  make  the  substitution  of  convenience  that  a  =  fo/T^.  Next,  we 
introduce  a  shifted  coordinate  system,  £  =  x  —  Sf,  which  allows  equation  5  to  remain  unaltered 
with  respect  to  £,  while  at  the  same  time  redefining  the  bounds  of  the  domain  as  0  <  £  <  A, 
where  A  —  Si  —  Sf. 

The  general  solution  to  equation  5  is 


u  =  Ae~°*  +  Beai  +  C  , 

(6) 

where  the  constants  A  and  B  are  subject  to  the  constraint 

2y/AB  =  -  . 

a 

(7) 

Note  that  equation  7  requires  that  A  and  B  be  of  the  same  sign,  so  that  their  product  is  positive. 
The  boundary  conditions  developed  in  section  3.3.1  require  that 

■u(O)  =  A  +  B  +  C  =  uq 

(8) 

u'(0)  =  (— A  +  B)a  =  —  tan  60 

(9) 

u(  A)  =  Ae~aA  +  BeaA  +  C  =  0 

(10) 

u\  A)  =  —Aae~aA  +  BaeoA  =  0  . 

(11) 
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Equation  11,  upon  multiplication  by  eaA  and  upon  elimination  of  a  by  way  of  equation  7,  may  be 
directly  solved  for  A  as 

A  =  VAB  In  (A/B)  .  (12) 

From  this  and  equation  7,  the  following  two  results  follow: 


e 


aA 


e 


— aA  _ 


i A/B 

4Wa 


(13) 


Substituting  these  results  into  equation  10  reveals  that 


C  =  -2y/AB 


(14) 


which  also  happens  to  equal  the  negative  reciprocal  of  a. 

Starting  with  equation  9  and  eliminating  a  with  a  substitution  of  equation  7  produces  the  equation 

yjA/B-  \[BjA  =  2  tan  0O  •  (15) 


This  may  be  solved  quadratically,  in  light  of  the  requirement  that  A  and  B  be  of  the  same  sign,  as 


1  +  sin  0O 
cos  do 


(16) 


We  are  now  prepared  to  substitute  equation  14  into  equation  8,  while  isolating  B,  to  obtain 


B  ((A/B)  +  1  -  2^A/B(j  =  uo 


(17) 


This  conveniently  reduces  to 

b(sJa/b-i}2  =  uo 


(18) 


The  term  B  may  now  be  directly  solved  by  substituting  equation  16  into  equation  18  in  order  to 
obtain 


B  =  M0(l  -  sin  dp) 
2(1  —  cos  6*o ) 


(19) 


The  quantity  A  directly  follows  as 


A  _  u0(l  +  sin  dp) 

2(1  —  cos0o) 


(20) 
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Now  that  A  and  B  are  explicitly  known  in  terms  of  the  boundary  conditions,  C  may  also  be 
explicitly  calculated  by  substitution  into  equation  14: 

(21) 


(22) 


These  parameters  may  be  employed  through  equation  6  to  trace  the  shape  of  the  cable  within  the 
process  zone,  where  the  distributed  force  /0  is  applied.  Figure  4  shows  the  actual  shape  of  the 
cable  for  different  load  angles  90.  For  smaller  load  angles,  the  curves  exit  the  right  side  of  the 
graph,  because  the  a:- measure  of  the  process-zone  size,  A,  exceeds  u0.  Therefore,  we  provide 
figure  5,  which  shows  the  typical  cable  geometry  for  various  load  angles,  except  that  the  x-axis 
has  been  normalized  by  A. 

Note  from  figure  5  that  the  normalized  shapes  of  the  curves  are  all  very  similar  for  load  angles 
below  60°.  As  the  load  becomes  more  perpendicularly  oriented  to  the  substrate,  however,  the 
shape  of  the  cable  acquires  much  more  of  a  knee-like  shape  to  it  and  the  tail  becomes  longer. 

4.2  Process-Zone  Cracklength 

Now  that  the  shape  of  the  deforming  surface  of  the  process  zone  has  been  fully  defined  through 
the  calculation  of  the  parameters  A,  B,  C,  and  A,  the  process-zone  cracklength  can  be  calculated. 
The  process  zone  is  defined  by  the  rectangle  Sf  <  x  <  Si  and  0  <  y  <  u0  (refer  to  figure  2).  It 
spans  an  x  distance  of  Si  —  Sf,  which  we  have  called  A  and  a  y  distance  of  u0.  Equation  22 
gives  us  an  explicit  measure  of  the  x-component  of  the  zone  size.  However,  A  is  not,  strictly 
speaking,  the  process-zone  cracklength.  From  figure  3,  we  note  that  the  length  of  the  crack 
formed  is  the  integral  of  ds,  not  that  of  dx.  Therefore,  the  full  process-zone  cracklength  s,  is 
obtained  as  follows: 

s  =  j  ds  =  J  Vl  +  u'2  dS,  .  (23) 

From  directly  differentiating  and  squaring  equation  6,  it  may  be  verified,  in  light  of  equation  7, 
that 

Vl  +  u'2  =  a(Ae~at  +  Beai)  .  (24) 


C  — 


-u  o 


sec  On 


The  quantity  A  follows  from  equations  12  as: 


A  = 


u  o 


2(sec0o  —  1) 


In 


1  +  sin  6*o ' 
1  —  sin  On 


Thus,  our  four  unknowns,  A,  B,  C,  and  A,  have  been  fully  characterized. 
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0.0  0.2  0.4  0.6  0.8  1.0 


5 /uo 

Figure  4.  The  process-zone  shape  (aspect  ratio  1:1),  with 
the  load  angle  do  as  a  parameter. 


0.0  0.2  0.4  0.6  0.8  1.0 


5 /A 

Figure  5.  The  process-zone  shape,  where  the  x-coordinatc 
has  been  normalized  by  A,  with  the  load  angle 
do  as  a  parameter. 
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Integrating  equation  24  across  the  process  zone  ( i.e .,  from  £  =  0  to  A),  in  light  of  equation  13, 
gives  the  following  value  for  the  cracklength  s  within  the  process  zone: 


A  comparison  of  s  to  A  is  presented  in  figure  6.  Naturally,  s  >  A  for  all  loading  angles  9q.  A 
specific  point  of  note  on  figure  6a  is  the  behavior  as  0o  approaches  90°  (i.e.,  applied  loading 
perpendicular  to  the  plane  of  symmetry).  We  note  that  the  a:- width  of  the  process  zone  (A) 
approaches  zero  length.  However,  the  pathlength  of  cable  in  the  process  zone  (i.e.,  the 
process-zone  cracklength  s)  can  never  fall  below  u0,  because  even  the  limiting  case  of  a 
process-zone  rectangle  approaching  dimensions  0  x  u0,  the  cable,  oriented  perpendicular  to  the 
plane  of  symmetry,  will  retain  a  cracklength  of  ?j0. 

Knowledge  of  A  allows  the  contents  of  figure  6  to  be  presented  in  a  slightly  different  manner. 
Namely,  tan  '(uq/A)  represents  the  average  angle  of  the  cable  within  the  process  zone.  This  is 
shown  in  figure  7. 

4.3  Bonding  Force  F0 

The  bonding  force  F0  is  also  of  interest.  We  define  it  as  the  total  resistive  force  in  the  ^/-direction 
applied  to  the  cable  in  the  process  zone  Sf  <  x  <  Si.  We  should  note  that,  because  of  the 
perpendicularity  of  F0  and  T^,  equilibrium  upon  the  cable  dictates  that  T02  =  F02  +  T^.  The 
resistive  force  comprising  F0  is  uniformly  distributed  along  the  cable  length  with  a  constant 
magnitude  of  /0  per  unit  cable  length.  Therefore,  calculating  the  total  bonding  force  is 
straightforward: 

Fq  =  fos  •  (27) 

Substituting  from  equation  26,  and  introducing  the  term  grouping 

O'  =  foUo  ,  (28) 

the  total  bonding  force  is  obtained  as 

F0  =  er  cot(#o/2)  .  (29) 
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(a) 


(b) 

Figure  6.  Normalized  x-component  (A/no)  and  pathwise 

(s/no)  size  of  process  zone  vs.  angle  of  applied  load 
6q,  expressed  on  a  (a)  linear  and  (b)  semi-logarithmic 
scale. 
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Figure  7.  The  averaged  cable  angle  within  the  process 
zone  vi.  the  load  angle  Oq. 


In  equation  28,  the  product  of  the  distributed  bonding  force  per  unit  cracklength  and  the  critical 
displacement  to  bring  about  full  fracture  represents  the  work  per  unit  cracklength  required  to 
completely  fracture  the  bond.  Thus,  a  is  a  material  property  that  may  be  termed  the  specific  bond 
energy  ( e.g .,  the  bond  energy  per  unit  cracklength). 

Therefore,  equation  29  shows  the  total  bonding  force  to  be  directly  proportional  to  the  specific 
bond  energy  a  and  the  cotangent  of  the  half-angle  of  the  applied  force.  In  the  limiting  case  of  a 
90  =  90°  applied  load,  the  total  bonding  force  (F0)  exactly  equals  the  specific  bond  energy  a. 

This  result  should  not  be  surprising  in  light  of  section  4.2,  which  showed  the  minimum  size  of  the 
process-zone  cracklength  to  be  u0  as  the  load  angle  90  approaches  90°.  Since  the  bonding  force 
F0  is  proportional  to  the  process-zone  cracklength  s  and  not  merely  the  x-component  of 
cracklength  A,  the  bonding  force  can  never  fall  below  the  specific  bond  energy  cr,  even  as  A 
approaches  zero  length. 

4.4  The  Relation  of  Applied  Tension  T0  to  Load  Angle  90 

The  bond  force  F0,  which  is  known  explicitly,  may  be  used  to  calculate  the  ultimately  desired 
quantity,  the  tension  in  the  cable  required  to  bring  about  separation  angle  90.  From  geometry,  we 
know  that,  in  order  to  satisfy  ^-equilibrium,  the  ^/-component  of  tension  must  exactly  equal  the 
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bonding  force  F0.  Thus, 


Therefore, 


Tn  = 


a 


sin0o  2  sin2  (6*0/2) 


sm  [2  r 


1  * 

2T0 


(30) 


(31) 


This  is  precisely  the  relation  obtained  in  reference  1  where  the  analogous  problem  was  solved 
using  variational  principles  of  energy  exchange  between  a  hanging  mass  m  and  the  specific 
chemical  energy  a  of  the  bond,  which  supports  the  mass  through  a  cable  bonded  to  a  rigid 
substrate.  In  that  solution,  the  tension  applied  to  the  bond  at  an  angle  90  arises  from  a  weight, 
given  by  mg.  However,  in  the  variational  approach,  the  mechanics  of  the  problem  are  not  even 
addressed  and,  thus,  it  was  a  bit  unexpected  (though  fully  welcomed)  to  achieve  exactly  the  same 
result. 

4.5  and  an  Interesting  Relationship 


The  quantity  T^  is  the  only  force  quantity  left  to  calculate.  It  is  easily  obtained  as 
^  =  T0  cos  6*(J,  which  reduces  to  the  following: 


T  = 

1  m 


cr 


sec  60  —  1 


(32) 


Figure  8  presents  the  relation  of  the  applied  cable  tension  T0,  the  total  bonding  force  T0,  and  T^ 
as  a  function  of  the  applied-load  angle  90.  Above  load  angles  of  45°,  F0  will  exceed  T^  in 
magnitude.  Below  45°,  the  reverse  is  true.  In  all  cases,  T0  exceeds  the  other  forces  because  it  is 
the  hypotenuse  associated  with  the  components  T0  and  T^. 

One  very  interesting  relation  that  may  be  shown  is  that 

T0  -  Tm  =  a  ,  (33) 

regardless  of  the  angle  of  applied  force,  60.  That  the  difference  of  T0  and  T^  should  remain 
invariant  may  be  understood  in  the  following  way.  Every  segment  of  cable  in  the  process  zone  is 
subjected  to  a  distributed  y-force  of  magnitude  /0.  For  each  segment,  there  is  a  component  of 
that  distributed  force  which  is  directed  along  the  cable  length,  thereby  changing  the  magnitude  of 
the  cable  tension.  This  distributed  component  directed  along  the  cable  is  of  magnitude  /0  sin  9, 
where  9  varies  from  0  to  90  through  the  process  zone.  The  overall  change  in  the  cable  tension, 
therefore,  is  /0S  /0  sin  9  ds.  The  quantity  /0,  being  constant,  may  be  taken  out  of  the  integral, 
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(a) 


(b) 

Figure  8.  Normalized  applied  force  T$/ <j  and  bonding  force 
Fo/cr  vs.  angle  of  applied  load  0 o,  expressed  on  a 
(a)  linear  and  (b)  semi-logarithmic  scale. 
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while  the  remaining  integrand  (sin  9  ds )  is  simply  du.  The  change  in  cable  tension  is,  therefore, 
/o  i(j'°  du  =  folio  =  a,  which  is  independent  of  the  particular  geometry  of  the  process  zone. 

Note  that  the  relationships  for  F0,  T0,  and  are  functions  of  the  term  grouping  a,  rather  than 
functions  of  f0  alone  or  u0  alone.  What  this  indicates  is  that,  while  f0  and  u0  were  introduced  as 
fundamental  model  parameters,  the  calculation  of  equilibrium  loads  is  really  not  dependent  on 
their  individual  respective  values,  but  rather  only  upon  the  load  angle  and  the  specific  bond 
energy,  a.  Where  u0  will  play  a  direct  role  is  in  the  calculation  of  the  field  energy  within  the 
process  zone  itself,  as  will  be  shown  in  section  4.6  of  this  report. 

4.6  Total  Field  Energy  of  the  Bond  and  the  Work  of  Deformation 

In  section  4.5,  the  specific  bond  energy,  a,  was  introduced  to  describe  the  energy  required  to 
create  a  unit  of  cracklength  (i.e.,  to  grow  the  crack).  We  also  have  the  ability  to  integrate  what 
would  traditionally  be  called  the  elastic  energy  of  the  bond.  In  the  current  model,  however,  the 
term  “elastic”  might  not  be  appropriate,  because  the  force-displacement  relationship  here  is  a 
constant,  and  not  linear  with  displacement  in  the  traditional  sense  of  elasticity.  Additionally,  the 
bonding  force  may  be  dissipative,  as  in  the  case  of  a  plastic  or  velcro  bond.  Therefore,  rather 
than  call  it  “elastic”  energy,  we  will  refer  to  it  as  the  total  field  energy  of  the  bond.  It  refers  to  the 
integration  of  energy  throughout  the  process  zone,  arising  from  the  displacement  of  the  cable 
from  the  u  =  0  configuration: 

E  —  j  fouds  =  /o  J  uVl  +  u'2  d£  .  (34) 

In  the  simplest  case,  if  the  cable  geometry  were  linear  within  the  process  zone  (rather  than 
catenoid),  such  that  the  cable  extended  linearly  from  the  (x,  y)  point  (Sj,  u0)  to  (St.  0),  the  total 
field  energy  of  the  bond  would  be  as/ 2. 

Except  when  the  load  is  applied  at  90°,  the  present  solution  does  not  have  a  linearly  varying 
displacement  in  the  process  zone.  For  the  special  case  of  90°  loading,  the  cable  approaches  in  the 
limit  a  y-oricntcd  line  in  the  process  zone,  of  cracklength  s  =  u0.  Let  us  define  a  baseline  field 
energy  as  e0  =  au0/2  =  j\]u//2.  It  represents  the  energy  stored  (or  dissipated)  in  the  process 
zone  of  a  bond  that  is  subjected  to  a  y-oricntcd  load  (i.e.,  perpendicular  to  the  plane  of  symmetry) 
of  limiting  magnitude  T0  =  F0  =  a. 

For  load  angles  other  than  90  =  90°,  we  would  expect  (because  of  the  displacement  field’s 
concavity)  the  total  field  energy  to  be  less  than  as/ 2.  Let  us  proceed  to  obtain  it.  Equation  34 
may  be  directly  integrated  through  substitution  of  equations  6  and  24,  in  light  of  equation  13. 
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The  result  is  given  by 


E 


fo 


A-B 


(. A  +  B  +  C)  + 


A-B 


c  +  aVab 


(35) 


which  may  be  reduced  to 

E  =  ^[as  -  T^s  -  A)}  ,  (36) 

consistent  with  our  expectations.  This  energy  relationship  is  displayed  in  figure  9,  normalized  by 

eo- 


To  understand  the  physical  origin  of  equation  36,  and  how  the  field  energy  arises  from  externally 
applied  work,  figure  10  is  provided.  The  point  O  represents  the  location  x  =  Si,  where  £  =  A. 
The  actual  path  representing  the  cable’s  deformed  state  is  given  by  the  arc  OB  (of  length  s).  If 
the  cable  lying  along  arc  OB  were  not  subject  to  the  tension  T0  at  point  B,  the  cable  would 
instead  lie  (unstressed)  along  the  line  OA. 


Figure  9.  Field  energy  (E/e o)  in  process-zone  bond 

vi.  load  angle  9q  (energy  as/2  normalized  by 
eo  given  for  reference). 


We  contemplate  the  application  of  tension  7^  to  point  A  of  the  unstressed  cable,  applied  at  a 
constant  angle  of  60  (in  this  case,  60°).  The  magnitude  of  the  tension  would  start  as  0  when  the 
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Figure  10.  Schematic  showing  the  origin  of  the  field 
energy,  in  terms  of  external  work  applied. 


cable  was  in  the  OA  configuration,  increasing  to  a  magnitude  of  T0  as  the  cable  reoriented  itself 
to  the  deformed  configuration  OB.  Because  the  magnitude  of  tension  T  can  be  shown  directly 
proportional  to  the  ^-displacement  of  the  endpoint,  one  may  assert  that  the  averaged  value  of 
tension  T  applied  when  the  point  A  is  deformed  to  point  B  is  exactly  (0  +  T0)/2  =  T0/2. 

The  work  W  done  by  this  action,  as  point  A  moves  to  B  under  increasing  force  ±  ,  is  the  integral 
of  the  force  times  the  component  of  displacement  in  the  direction  of  action,  which  is  /  d AB, 
given  more  simply  as  W  =  (T0/2)\CB\.  Because  OE  is  parallel  with  7\  while  AE  and  GF  are 
perpendicular  to  it,  we  can  assert  that  \CB\  =  \DB\  —  \OE\  +  | OF\.  We  may  further  assert,  on 
geometrical  grounds,  that  LFOG  =  Z BGD  =  90.  Therefore,  \DB\  =  \GB\  sin 0O  =  u0sin90. 
Likewise,  | OE\  =  \OA\  cos  60  =  s  cos  90-  Finally,  we  have  that  \OF\  =  \OG\  cos  60  =  A  cos  60. 
Substituting  all  these  metrics  into  the  equation  representing  the  work  done  to  deform  the  cable 
from  the  unstressed  configuration  OA  to  the  fully-stressed  configuration  OB  gives  the  following 
result: 

W  =  (T0/2)[m0  sin  90  —  (s  —  A)  cos  6*0]  .  (37) 
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Knowing  that  T0u0  sin  90  =  F0u0  =  as  and  that  T0  cos  90  =  allows  equation  37  to  be  reduced 
to 

W  =  ^[as-T00(s- A)]  .  (38) 

The  identical  forms  of  equations  36  and  38  tell  us  that  the  amount  of  field  energy  stored  in  (or 
dissipated  by)  the  bonds  along  the  cable  in  the  region  OB  exactly  equals  the  external  work  that 
would  have  been  required  to  deform  the  cable  from  its  stress-free  configuration  OA  to  its  fully 
loaded  configuration  OB. 

4.7  Axial  Cable  Displacement  with  Changes  in  90 

Examining  the  underlying  premise  of  figure  1,  it  is  apparent  that,  as  the  cable  bond  becomes 
progressively  “unzipped,”  the  weights  (each  of  adjustable  mass  m)  will  give  up  potential  energy 
as  they  lower  in  height.  In  this  section,  we  seek  to  determine  the  amount  of  axial  cable 
displacement  (i.e..  decrease  in  height  of  mass  m )  that  occurs  as  a  function  of  the  load  angle  90. 

As  a  reference  state,  we  chose  the  situation  where  the  initial  value  of  90  is  90°  and  the  initial  cable 
tension  exactly  equals  the  associated  equilibrium  value  of  T0  =  a.  As  the  masses  lower  and  the 
cables  begin  to  debond,  the  decrease  in  cable  length  between  points  H  and  P  corresponds  exactly 
to  the  change  in  the  elevation,  Az,  of  the  masses. 

We  employ  figure  1 1  to  assist  in  determining  cable  lengths.  Initially,  we  assume  that  u0  (and  thus 
the  process  zone)  is  vanishingly  small,  so  that  the  fractured  geometry  associated  with  load  angle 
90  is  composed  of  two  segments:  the  detached  cable  segment  (of  length  c)  and  the  bonded  cable 
segment  (of  length  p  —  a).  Under  this  assumption,  the  elevation  change  of  the  masses,  Az,  when 
the  load  angle,  90,  deviates  from  90°,  is  given  by 

Az  =  h  +  a  —  c  (when  no  h)  ,  (39) 

where  a  =  h/  tan  9{i  and  c  =  h/  sin  9(i. 

However,  when  the  size  of  the  process  zone  is  not  considered  negligible,  a  further  refinement  to 
equation  39  is  required  to  obtain  the  elevation  change  of  the  masses,  Az:,  when  the  load  angle,  90, 
deviates  from  90°.  Employing  the  relation  for  s  from  equation  26  and  seeing  that  b  =  uq/ sin  90 
and  d  =  A  —  n0/  tan  90,  we  obtain 

h 

Az  =  (h  +  a  —  c)  +  (b  +  d  —  s)  =  - —  (sin  90  +  cos  90  —  1)  +  (A  —  2 u0  cot  90)  .  (40) 

sin  90 
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Figure  1 1 .  Schematic  showing  lengths  associated  with  cable  displacement. 


To  repeat,  the  value  of  Az  represents  the  finite  change  in  elevation  level  of  the  masses  of  figure  1 
as  a  function  of  the  load  angle  60.  An  initial  value  of  90  =  90°  is  assumed.  The  first  term  of 
equation  40,  proportional  to  h,  is  the  dominant  term,  arising  from  the  macroscopic  geometry  of 
the  system.  The  second  term,  on  the  order  of  u0,  provides  a  secondary  correction,  to  account  for 
the  catenoid  shape  of  the  process  zone.  Further  insight  could  be  pursued,  for  example,  by 
comparing  the  displacement,  A z,  of  the  adjustable  masses  with  the  debonded  cable  length,  given 
by  Sf  =  ( h  —  Mo)/ tan  0o- 

4.8  The  Question  of  Asymptotes 

One  of  the  questions  that  arose  (in  the  minds  of  the  authors)  during  the  course  of  this  investigation 
was  whether  or  not  the  eventual  solution  would  be  asymptotic.  That  is  to  say,  whether  the 
location  of  St  (see  figure  2)  would  approach  infinity,  even  for  finite  values  of  T0.  The  present 
solution  shows  that,  for  the  current  set  of  governing  assumptions,  it  does  not.  One  mathematical 
fact  regarding  asymptotes  is  that  their  second  derivative  approaches  zero,  since  an  asymptote’s 
graph  approaches  a  straight  line. 

Considering  equation  5,  however,  one  can  see  that,  as  long  as  the  specific  bonding  force  f0  is 
nontrivial,  u"  must  be  greater  than  zero.  Thus,  the  governing  equation  5  is  incapable  of 
producing  an  asymptotic  result,  wherein  u"  =  0. 


20 


This  also  raises  the  possibility,  however,  that  were  the  bonding-force  relation  something  other 
than  a  constant,  the  conclusion  regarding  the  possibility  of  an  asymptotic  solution  would  change 
as  well.  For  example,  instead  of  a  constant  specific  bond  force  /0,  the  bond  might  be 
characterized  by  an  spring-type  f  =  ku  specific  bonding  force.  With  such  a  physical  model,  the 
governing  relation  would  become  u"  —  {k/T^u^l  +  u'2  =  0.  In  this  circumstance,  it  would  be 
mathematically  possible  to  achieve  an  asymptotic  solution,  in  which  u  approached  zero  only  as  x 
approached  infinity.  In  such  a  case,  u"  would  likewise  approach  zero,  which  is  the  requirement 
for  an  asymptote. 


5.  Conclusions 


In  this  report,  a  notional  experiment  was  presented  involving  the  progressive  debonding  of  two 
flexible  surfaces  (e.g.,  tape  or  velcro).  An  idealized  model  was  developed  which,  when  solved, 
predicts  the  shape  of  the  flexible  substrate  in  the  process  zone  where  debonding  occurs.  More 
importantly,  it  provides  the  force  and  energy  relationships  that  govern  the  problem  as  a  function  of 
the  loading  angle,  which  will  provide  the  means  to  test  the  model  against  future  experimental  data. 

The  assumptions  contained  in  the  formulation  are  that  (1)  the  (tape  or  velcro)  substrate  is  fully 
flexible;  (2)  the  bonds  are  one-dimensional  in  nature,  such  that  there  is  no  lateral  (;. e.,  Poisson) 
interaction  in  the  process  zone;  (3)  the  bonding  force  in  the  process  zone  is  independent  of 
separation  up  to  the  point  of  fracture  (in  the  manner  of  a  perfectly -plastic  bond);  and  (4)  the  bonds 
are  uniformly  distributed  along  the  substrate  length. 

Based  on  these  assumptions  to  the  problem,  the  process  zone  was  seen  to  take  on  the  shape  of  a 
catenary  (i.e.,  a  free-hanging  cable  under  the  influence  of  gravity),  for  which  an  analytical 
solution  presents  itself.  Importantly,  the  force  relations  predicted  for  this  model  identically 
matched  those  developed  for  the  problem  using  an  alternate  energy-based  approach  (7). 
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Appendix.  The  Catenary  Solution 


The  catenary  is  a  textbook  problem  in  the  teaching  of  differential  equations.  Historically,  it  refers 
to  the  shape  a  chain  (or  cable)  assumes  under  the  influence  of  gravity,  being  attached  only  at  its 
endpoints.  We  rederive  the  result  here  for  the  convenience  of  the  reader. 

The  governing  equation  to  solve  is  functionally  identical  to  that  given  in  equation  5,  namely: 

u  —  aV  1  +  u'2  =  0  .  (A-l) 


We  reduce  the  equation  order  by  casting  u"  as  dv! /du  ■  du/dx ,  knowing  that  du/dx  is  the  same  as 
v! .  In  so  doing,  we  obtain 


'1  +  ur 


=  adu  . 


This  is  readily  integrated  as 


'l  +  u'2  =  au  +  Ci  . 


(A-2) 


(A- 3) 


Now  we  note  that,  in  equation  A-l,  a  vertical  shift  of  magnitude  wshift,  such  that  u  =  u  —  wshift 
produces  the  identical  equation  in  u.  This  indicates  that  any  solution,  u(x)  is  likewise  a  solution 
u(x).  The  vertical  shift  is  related  to  constant  C\.  Therefore,  we  may  choose  any  particular  value 
of  Ci,  for  mathematical  convenience,  knowing  that  any  solution  subsequently  obtained  for  u(x) 
can  be  subjected  to  an  arbitrary  vertical  shift.  We  choose  C i  so  that  u'  is  identically  zero,  when  u 
is  zero.  Such  a  choice  forces  C\  to  a  value  of  1.  With  this  substitution,  equation  A-3  becomes 

- ,  =  adx  .  (A-4) 

\fuJu  +  2/ a 


Introducing  the  substitution  v2  =  u  reduces  the  subsequent  integral  of  equation  A-4  to  a 
convenient  form: 


(A- 5) 


23 


The  left-side  integral  may  be  solved  in  the  following  way: 


f  dx 

f  (x  +  \/  x2  +  b2)dx  j 

1  H — ,  dx 

'  V  Vx2  +  b2J 

1  \/x2  +  b  J 

(x  +  \/x2  +  b2)Vx2  +  b  J 

x  +  \/  x2  +  b2 

=  ln(x  +  \Jx2  +  b2)  +  C 

With  this  integration  form,  equation  A- 5  may  be  integrated  as 

ln(n  +  yjv2  +  2/a)  =  (a/ 2)x  +  in  C2 

Equation  A-7  may  be  explicitly  solved  for  v  to  yield 

_  ^_pax/ 2  _  _J_  -ax/2 

2  aC2 

Squaring  this  equation  and  substituting  u  for  v 2  gives  the  following: 

u=T^  +  Jcl^-1/a  ' 

Equation  A-9  is  equivalent  to  the  form 


(A-6) 


(A-7) 


(A-8) 


(A-9) 


u  =  Ae~ax  +  Beax  +  C  ,  (A- 10) 

where  A  and  B  are  subject  to  the  constraint  2  \/A  B  =  1/a.  We  see  that  the  solution,  given  by 
equation  A- 10,  matches  the  form  presented  in  the  report  as  equation  6,  subject  to  the  identical 
constraint  on  A  and  B.  Often,  in  the  literature,  the  result  is  given  in  terms  of  the  equivalent 
hyperbolic  cosine  function. 

For  the  general  catenary  solution  (which  may  be  shifted  vertically),  C  is  a  free  constant.  If  the 
constraint  is  enforced  to  make  v!  =  0  when  u  =  0  (as  was  enforced  when  we  chose  6)),  then  C 
will  take  on  a  value  of  —1/a,  as  dictated  by  equation  A-9.  This  constrained  value  of  C  was 
indeed  obtained  in  the  report  as  equation  14,  in  which  C  was  revealed  as  —  1/a. 
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